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We report the results of calculation of diffusion coefficients for Lennard-Jones 
liquids, based on the idea of time-scale invariance of relaxation processes in liquids. 
The results were compared with the molecular dynamics data for Lennard-Jones 
system and a good agreement of our theory with these data over a wide range of 
densities and temperatures was obtained. By calculations of the non-Markovity 
parameter we have estimated numerically statistical memory effects of diffusion in 
detail. 
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The motion of atoms and molecules in liquids is a typical example of a many-body prob- 
lem. It is well known, however, that in contrast to gases and solids the physical mechanisms 
of transport processes in liquid matter is no well established In particular, for diffusive 
systems the role of memory effects induced by the disorder of the background medium is not 
investigated enough There are different theoretical approaches and approximations 

for calculating transport coefficients JUjU | , in some of them they resort to time correlation 
functions (see, for example, I). The general approach of time correlation func- 

tions leads to calculation of exact expressions of transport coefficients for hydrodynamic 
equations. The similar expressions are known as the Green-Kubo relations. 



In the present work we suggest a new approach to calculating transport properties 
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on the realization of idea of time-scale invariance of relaxation processes in liquids 3| by 
well-known Zwanzig- Mori's memory function formalism 0,0]. Recently this idea has made 
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it possible to explain the experimental data on slow neutron scattering in liquid cesium 
and sodium 0]. Here we have tested this approach on Lennard- Jones (LJ) liquids. As 
known, among a large set of dense fluids these systems are one of the most popular (base) 
for different investigations. So, diffusion coefficients of LJ liquids were calculated for a wide 
range of densities and temperatures, and then quantitative analysis of the memory effects in 
diffusion phenomena has executed. The obtained results are compared with the predictions 

n n 

of other theories j6| and the molecular dynamics data |10| . 

As known, the diffusion coefficient can be expressed in terms of the velocity autocorrela- 
tion function (VACF) 



by the relation 



hr^T roc 

D = ^ a{t)dt. (2) 
m Jo 

The latter is one of the Green-Kubo equations, relating the integrals of the correlation 

n 

function to the transport coefficients Here kB, T and m are the Boltzmann constant, 
temperature and atomic mass, respectively. 

On the other side, the generalized Langevin equation (GLE) for VACF can be written in 
the following way: 

'^''^^^ -uj^^U' M,{r)a{t-r)dT, (3) 



dt JO 

where Mi(t) is the first-order normalized memory function [Mi(t = 0) = 1] and cu^"^^ is the 
second frequency moment of VACF. This equation is obtained by the use of Zwanzig- Mori's 
projection operators formalism However, the projection operators technique allows to 
derive the whole chain of interconnected equations similar to the above mentioned equation. 
These equations contain the memory functions Mj(t) and the even frequency moments u;^^*^ 
of higher orders, i.e. i = 1,2,3,.... This chain can be represented in terms of Laplace 
transforms as a continued fraction: 



ais 



dte-''a{t) = [s + cu(2)Mi(s)]'1 = V + cu^^)/ [s + (c^(^V^^'^ - c^^'^)M2(s)] } 
1/ {s + ^(2)/ [s + - + •••)] } , (4) 
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and Eq. Q yields 

D = ^d{s = 0). (5) 
m 

So, the problem of calculation of transport properties can be reduced to calculation 
of the certain memory functions Mi{t). Although there are microscopic expressions and 
some formal prescriptions for memory functions calculations, it is prohibitively difficult for 
models like Lennard- Jones (LJ) liquid. A rich variety of model functions was proposed for 
this purpose, most of which have little physical justification. For example, hyperbolic secant 
memory was used to calculate velocity, transverse stress , energy current density correlation 



functions and the correspond.ng transport coefficients for LJ hqn.d fl. The val.dity and 

justification of this approach was closely examined in work 

There is another more effective and powerful way to calculate transport coefficients of 



121 to calculate self- 



the diffusion constant type. This method was first suggested in work 
diffusion in liquid argon, and used later in Ref. [l^ to investigate viscosity effects in liquid 
argon. Here we use this approach to analyze diffusion time-scales and to calculate the 
diffusion coefficient in LJ liquid over a wide density and temperature ranges. 

As known, the memory functions Mj(t) have characteristic time scales, which can gener- 
ally be defined at fixed z by = Mj(s = 0) = /q°° dtMi{t), where Tj is the relaxation time of 
a certain correlation function Mi{t). These time scales characterize the corresponding relax- 
ation processes and can have different numerical values. Nonetheless, on a certain relaxation 
level, (for example, on the ith level) the scale invariance of the nearest interconnected re- 
laxation processes can exist. We shall consider below the following approximation Xj+i = Xj. 
Physically it means the occurrence of time-scale invariance on the ith relaxation level. 

In case of i = we have an approximate equality of VACF relaxation time Tq and the 
relaxation time of the first-order memory function ti. Then, from Eq. (j^ and the second 
equality of Eq. (jlj) at ti = tq one can obtain 



m 



(2) 



1/2 47rfi 

, ur' = / drg{r)r'^ 

3 Jq 



3dU(r) d (dUir) 



r dr dr \ rdr 



(6) 



where n is density, g{r) is radial distribution function and U{r) is interparticle potential. 

At i = 1 we obtain the equality of the relaxation times of first and second order memory 
functions In this case, T2 = ti, and from Eqs. Q and ^ we have 
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Figure 1: The reduced density n* = na^ dependence of diffusion coefficient at the reduced tem- 
perature T* = kBT/e = 3.46, 2.5, 1.81 and 1.23. The triangles (A A A) show results of Eq. ®, 
the solid line presents our results from Eq. the circles (o o o) show results of Eq. ©, the 
dotted line (• • •) presents the results of approximation with hyperbolic secant memory of Ref. Q], 

the broken ( ) and the chain line (— • — ) correspond to results of two models for diffusion 

coefficient obtained on the basis of Joslin and Gray study in terms of the first two and three Mori's 
coefficients, correspondingly P, The full circles (•) represent molecular dynamics data of Ref. 
|lol | and correspond to diffusion values in points well within the fluid region of the LJ phase dia- 
gram 16]. Presented data bring out clearly that results of our theory are in excellent agreement 
with molecular dynamics data without any deviations. 
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where g3{r, ri) is the triplet correlation function, and f3 is the cosine of the angle between r 
and Ti. 

Following the same procedure for third and second order relaxation times, T3 = T2 (and 
i = 2), we come to the following expression: 

m [cj{2)]3/2[c^{6)^(2) _ (^(4))2]l/2- W 

In a more general case with = rj_i, i = 1, 2, ... we can obtain the expression of diffusion 
coefficient by the first even frequency moments: uj^^\... ,uj^'^'^\ The second, the forth 

and the sixth frequency moments contained in Eqs. © and (jHI) were approximately 

obtained for LJ liquid by authors [6]. Using these values, the diffusion coefficients were 
calculated with the help of Eqs. (jH)), (0) and (jSI). The comparison shows that the values of 
D obtained by Eq. ((Zj) have the best agreement with the molecular dynamics data for the 
whole studied range of densities and temperatures. The results for D* = D{m/ ea'^Y^'^ from 
Eq. (jH)) are presented in Fig. 1 as triangles, from Eq. (|7j) are shown as solid curves^ and 
from Eq. (jSJ are presented by circles. We also represented here the results of authors [6|] for 
comparison with our ones. Namely, the diffusion coefficient obtained by Ref. 0] from the 
approximation of the second-order memory function by the hyperbolic secant are shown as 
dotted curves, whereas those obtained in Ref. from the prescription of Joslin and Gray 
jisl ] are shown by broken and chain curves. From Fig. 1 one can easily see that the results of 
Eq. ((7j) have a better agreement with the molecular dynamics data jiol . Q| presented by full 
circles in more cases than data found of others theories. This result is of great interest for 
the study of diffusion phenomena. It confirms the possibility of equality of relaxation times 
of first and second order memory functions, i. e. Ti and T2. As can be seen from Fig. 1, the 
approximation with hyperbolic secant memory show a sufficiently satisfactory agreement 
for some cases. Joslin and Gray models yield an understated results in comparison with 
molecular dynamics data. The distinction between these theories and molecular dynamics 
data grows with decreasing of density. However, all examined models, including even Eqs. 
(jni) and (jHI), begin to reproduce qualitatively molecular dynamics data as the triple point of 
LJ system with the reduced parameters n* = 0.849, T* = 0.773 is approached. 

The next step of our study is estimation of statistical memory effects in diffusion pro- 
cesses of LJ liquid. It can be successfully done within the framework of memory function 

I I 

and dimensionless non-Markovity parameter formalism. The last one is the ratio of 
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Figure 2: Variation of the non-Markovity parameter eq for VACF of LJ liquid with temperature 
at different densities. Correspondence between the curves and the densities is presented in the 
inset. The behavior of the parameter eq indicate a quantitative change of ratio between VACF 
time scale and memory relaxation time. The cross-over from Markovian (eq ^ 1) to a quasi- 
Markovian (eq > 1) character of the diffusion process at low temperature and high densities is 
noteworthy. Referring to the obtained data the diffusion process is in fact quasi-Markovian in a 
wide temperature region at high densities. 

the relaxation time Tq of initial correlation function and the relaxation time Ti of memory 
function 

^0 = ro/n. (9) 
Using Eq. for the relaxation time Tq = a{s = 0) we obtain 

To 

From the second equality in Eq. (^j) at g 
expression 



and Eqs. © and fTIH) we come to the final 
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Figure 3: Variation of the non-Markovity parameter eq with the reduced density at four different 
temperatures T* = 1.23, 1.81, 2.5 and 4.5 reveal a nonlinear character of demarkovization at matter 
densifying. The horizontal dash line in the bottom of the figure corresponds to the quantitative 
value e = 1 (a full non-Markovian relaxation scenario). The markovization of the process occurs at 
decrease of density and increase of temperature; it is accompanied by increase of the parameter eq 
and partial disordering of LJ system. Increase of density results in ordering and demarkovization 
of diffusion process with sharp reduction of eo values. 
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On the basis of the definition © it is possible to define the spectrum of non-Markovity 
parameter: Ei = Ti/ri^^i (see Ref. [3l for a more detail). Then the general form of this 
parameter for i = 1,2, 3, ... can be expressed in terms of Mori's coefficients (static correlation 
functions) as 



'0 ' 



'+1,2 



if i is odd, 
if i is even. 



(12) 



Obviously, the parameter ei obtained by Eq. ()12p estimates the memory effects of the 
ith level relaxation process, which is described by the memory function Mj(t). 
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To investigate the memory effects of VACF itself it is necessary to study the zero point in 
spectra of the non-Markovity parameter Eq calculated by Eq. (fTT|) . Therefore, the parameter 
Eq has been calculated for a wide range of densities and temperatures for LJ liquid. The 
results of our calculations are presented in Fig. 2. From this figure one can see that the 
non-Markovity parameter always satisfies condition Eq > 1. This is the evidence of weak 
statistical memory effects in diffusion processes, which is usually observed in Markovian and 
quasi-Markovian processes. In this case ti is much smaller than the VACF relaxation time 
tq. The received result helps to understand the reason of amazing efficiency of Markovian 
approximation in the analysis of diffusion in LJ liquids (see, for example, one of the recent 
investigations ^^). It also enables to estimate the source of unsatisfactory agreement of 
Eq. (jni) with molecular dynamic data. The point is that Eq. ^ is based on the assumption 
of Ti = Tq and/or Eq = 1, which is true for non-Markovian processes only, while Eqs. (|7j) 
and (jHl) are applied both to Markovian and non-Markovian phenomena. Furthermore, in 
Fig. 2 we can see, that the parameter Eq smoothly increases with the increase of the reduced 
temperature T* = ksT/E, at the same time it decreases with the increase of density. So, 
non-Markovian effects fall off (and/or Markovity is enhanced) when the temperature of the 
system increases while the value of the density n* = na^ is fixed, i. e. markovization 
of processes is observed, and it reflects the amplification of randomness during diffusion. 
However, a "saturation " is attained at a certain value of temperature, and the values of 
the parameter Eq almost cease to change in this case. As may be seen from Fig. 2, such 
"saturation" is observed at lower temperatures for a more dense medium. The density n* 
dependence of the parameter Eq is shown in Fig. 3 for the four isotherms. The behavior of 
these curves has a nonlinear character. In this figure we can observe demarkovization, that 
is, amplification of regularity and robustness in the system which appears with the increase 
of density. For example, for the isotherm T* = 1.81 the parameter Eq decreases more than 
five times within density range 0.2 ^ 0.7. Non-Markovian effects are enhanced and begin 
to dominate in the vicinity of the triple point of LJ liquid, where a well-known negative 
correlation in the behavior of VACF is observed. The ratio between VACF time-scale and 
relaxation time of memory in this phase region is ~ 2.5 ^ 5. This is a quantitative evidence 
of considerable memory effects incipient diffusion process of LJ system near the triple point, 
where the density variation is short ranged and of order of a couple of molecular diameters 

Q 
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The results of this work can be summarized as follows, (i) We have presented a new 
approach in calculation of transport coefficients. The approach is based on the time-scale 
invariance idea, developed within framework of Zwanzig-Mori's formalism. Our theory allows 
to calculate transport properties in terms of frequency moments of corresponding correlation 
functions without any adjustable parameters, (ii) Three different expressions were obtained 
for the diffusion coefficient of LJ liquid and tested together with other theories and the 
molecular dynamics data over a wide range of densities and temperatures. The comparison 
showed that the equation obtained from the condition of approximate equality of relaxation 
times of the first- and second-order memory functions has the best agreement with the 
molecular dynamic data. This equation includes only the second and the fourth frequency 
parameters, which were calculated with a high degree of precision, (iii) The non-Markovity 
parameter calculated for the diffusion phenomena in LJ liquid reveals a Markovian character 
of thermal motions of particles. Finally, the values of this parameter demonstrate the cross- 
over from Markovian to quasi-Markovian relaxation scenario at low temperatures and high 
densities and allow to estimate quantitatively this cross-over. 
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